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^ I 1 INTRODUCTION 

^ ' While only a handful of directly-detected exoplanets is cur- 
rently known, the near future should bring a statistically sig- 
nificant sample of directly-imaged exoplanets, thanks to a 
number of surveys underway or coming online soon. Ex- 
amples include VLT/SP HERE , Gemini/GPI, Subaru/HiCIAO, 
Hale/P1640 j vigan et al.| llOlOl: Ich auvin et al. 2010; Liu et al. 



ABSTRACT 

The initial entropy 5i of a gas giant planet is a key witness to its formation history and a 
crucial quantity for its early evolution. However, formation models are not yet able to reliably 
predict Si, making an observational determination of initial entropies essential in order to 
guide formation scenarios. Using a grid of models in mass and entropy, we show how to place 
joint constraints on the mass and initial entropy of an object from its observed luminosity 
and age. Moreover, we demonstrate that with mass information, e.g. from dynamical stability 
analyses or radial velocity, tighter bounds can be set on the initial entropy. We apply this 
procedure to 2M1207 b and find that its initial entropy is at least 9.2 A;B/baryon, assuming that 
it does not burn deuterium. For the planets of the HR 8799 system, we infer that they must have 
formed with Si > 9.2 fce/baryon, independent of the age uncertainties for the star. Finally, a 
similar analysis for /3 Pic b reveals that it must have formed with Si > 10.5 fce/baryon, using 
the radial-velocity mass upper limit. These initial entropy values are respectively ca. 0.7, 0.5, 
and 1.5 ^a/baryon higher than the ones obtained from core accretion models by Marley et 
al., thereby quantitatively ruling out the coldest starts for these objects and constraining warm 
starts, especially for /3 Pic b. 

Key words: planets and satellites: general - methods: analytical - techniques: imaging spec- 
troscopy - stars: individual: 2MASSWJ 1207334-393254- stars: individual: HR 8799 - stars: 
individual: /3 Pictoris. 
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2010"; McBride et a]J 120111 : iMcElwain et all |2012| : iHinklev et alj 
201 lb) An important input for predicting the yields of such sur- 



veys and for interpreting their results is the cooling hi s tory o f 
gas giant plan e ts. In t he traditional approa ch Istevensonl ( Il982h : 
iBurrows et al] ( Il997l) : ISaraffe et al] j2003l) . objects begin their 
cooling, fully formed, with an arbitrarily high (specific) entropy 
and hence radius and luminosity. In the past, the precise choice of 
initial conditions has been of no consequence because only the 4.5- 
Gyr-old Solar System's gas giants were known, while the high ini- 
tial entrop y is forg ott en on the short Kelvin- Helmholtz timescale 
GM-/RL ( IStevensoi]ri982l : iMarlev et'alj|2007i) . However, direct- 
detection surveys are aiming specifically at young (^ 500 Myr) 
systems, and the traditional models, as their authors explicitly 
recognised, are not reliable at early ages. 

Using the standard core-accretion formation model 
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jPollack et al.1 Il99d : iBodenheimer et alj l2000l : iHubickvi etal] 
'2005*). M arley et al.l te007'. hereafter iMOVh found that newly- 
formed gas giants produced by core accretion should be substan- 
tially colder than the usual cooling tracks that begin with arbitrarily 
hot initial conditions assume. Thes e outc omes are known as 'cold 
start' and 'hot start', respectively. ImotI however noted that there 
are a number of assumptions and approximations that go into the 
core accretion models that make the predicted entropy uncertain. 
In particular, the accretion shock at the surface of the planet is 
suspected to play a key role a s most of the mass is processed 
through it (M07; Baraffe, Chab rier. & B arman 2010) but there 
does not yet exist a satisfactory treatment of it. Furthermore, 
there may also be an accretion shock in the other main formation 
scenario, gravitational instability, such that it too c ould yield cooler 
planets than usually expected (see section 8 of iMordasini et al] 
.2012b ). 

The most reasonable viewpoint for now is therefore to con- 
sider that the initial entropy is highly uncertain and may lie almost 
any where between the cold values of M07 and the hot starts. In 
fact, ImOTI calculated 'warm start' models tha t were intermediate 
betwe en the cold and hot starts, and recently, ISpiegel & Burrow^ 
( l2012h calculated cooling tracks and spectra of giant planets begin- 
ning with a range of initial entropies. This uncertainty in the initial 
entropy means that observations are in a privileged position to in- 
form models of planet formation, for which the initial entropy (or 
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luminosity) is an easily accessible output. Its determination both on 
an individual basis as well as for a statistically meaningful sample 
of planets would be very valuable, with the latter allowing quanti- 
tative compariso ns to planet population synthesis jlda & LinI 20041 : 
iMordasini et al . 2012b a). Contrasting with the current small num- 
ber of direct detections, which almost exclusively consist of mas- 
sive objects on wide orbits, the near future should bring both qual- 
itative and quantitative changes to the collection of known ones. 
Indeed, upcoming and ongoing surveys such as SPHERE, GPI, Hi- 
CIAO, JWST are expected to reveal a number of light objects at 
small separations from their parent star (i.e. Jupiter-like planets) 
It is therefore crucial to highlight now the need to interpret direct 
observations in a model-independent fashion and to think about in- 
formation these yield about the initial conditions of gas giants. 

In this paper, we investigate the constraints on the mass M and 
initial entropy S\ that come from a luminosity and age determina- 
tion of a directly-detected object, focussing on exoplanets. Since 
the current luminosity increases both with M and with S\, these 
constraints take the form of a 'trade-off curve' between M and 
Si. We show that, in the planetary regime, the allowed values of 
M and 5; can be divided into two regions. The first is a hot-start 
region where the initial entropy can be arbitrarily high but where 
the mass is well determined. This corresponds to the usual mass 
determination by fitting hot-start cooling curves. The second cor- 
responds to solutions with a lower, tightly-constrained initial en- 
tropy for which the planet mass has to be larger than the hot-start 
mass. For masses above the deuterium-burning limit, the required 
entropies drop more steeply with increasing mass due to energy 
injection by deuterium burning, but we defer a discussion of this 
to a forthcoming paper. The degeneracy between mass and initial 
entropy means that in general the mass and entropy cannot be con- 
strained independently from a measurement of luminosity alone. 
When additional mass information are available, for example for 
multiple-planet systems or systems with radial velocity, it is possi- 
ble to constrain the formation entropy more tightly. 

We start in Section|2]by describing our gas-giant cooling mod- 
els, discussing the luminosity scalings with mass and entropy, and 
comparing to previous calculations in the literature. In Section [3] 
we show in general how to derive constraints on the mass and ini- 
tial entropy of a directly-detected exoplanet by comparing its mea- 
sured luminosity and age to cooling curves with a range of initial 
entropies. We briefly consider in Section|4]similar constraints based 
on a (spectral) determination of the effective temperature and sur- 
face gravity. We apply the procedure based on bolometric luminosi- 
ties in Section|5]to three particularly interesting systems - 2M1207, 
HR 8799, and j8 Pictoris - and derive lower bounds on their initial 
entropies. For HR 8799 and /? Pic, we use the available information 
on the masses to constrain the initial entropies more tightly. Finally, 
we offer a summary and concluding remarks in Section |6] 



2 COOLING MODELS WITH ARBITRARY INITIAL 
ENTROPY 

Given its importance in the initial stages of a planet's cooling, 
we focus in this paper on the internal entropy of a gas giant as a 
fundamental parameter which, along with its mass, determines its 
structure and controls its evolution. Very few authors have explic- 
itly considered entropy as the second fundamental quantity even 
though this approach is more transparent than the usual, more in- 
tuitive use of time as the second independent variable. In this sec- 
tion, before comparing our models to standard ones, we discuss in- 



terior temperature-pressure profiles and semi-analytical ly explain 
the different scalings of luminosity on mass and entropy which ap- 
pear in the models. 

2.1 Calculation of gas giant models 

We calculate the evolut ion of co oling gas giants by a 'following 
the adiabats' ap proach jHubbard 1977; Fortney & Hubbard 2004; 
lArras & Bildsten 2006). We generate a large grid of planet models 
with ranges of interior specific entropy S and mass M and deter- 
mine for each model the luminosity L = L(M, 5 ) at the top of the 
convection zone. We then use it to calculate the rate of change of 
the convective zone's entropy AS lAt, which is defined by writing 
the entropy equation dLjdm = -TdS jdt as 

d5 r 

L = rdm + Lo, (1) 

dt J 

with Ld the luminosity due to deuterium burning in the convection 
zone. The assumptions in this approach are that dS /dt is constant 
throughout the convective zone and that no luminosity is generated 
in the atmosphere. For these t o hold, we require both the thermal 
timescale ftheira('') = PcpT/gF ( I Arras & Bildstej 2006) - where Cp 
is the specific heat capacity, g the gravity, and F the local flux - in 
the radiative zone overlying the convective core and the convective 
tumovef] timescale to be much shorter than the timescale on which 
the entropy is changing, ts = -5/(d5/df) = MTS/L, with f the 
mass-weighted mean temperature in the convection zone. A simi- 
lar approach was used to study the evolution of o hmically-heated 
irradiated gas giants bv lHuang & CummingI (l2012h . 

We do not include deuterium burning directly in the planet 
models since it always occurs inside the convection zone. Instead, 
we calculate the deuterium burning luminosity per unit deuterium 
mass fraction for each model in the grid and use it in equation (TJ to 
follow the cooling of the planet and the time evolution of the deu- 
terium mass fraction averaged over the convective region. We use 
the deuterium burning rate from Caughlan & Fowler (1988|) and 
includ e the inter mediate screening co r rectio n from lOraboske et alj 
Il973*), following [Burrows & Lieberll (Il993h . We adopt a standard 
initial D/H number ratio of 2 x 10"' (i.e. a deuterium mass fraction 
Xo = 3.0 X 10"'). Since we focus here on directly-detected planets, 
we defer a detailed exploration of deuterium burning in our models 
to an ulterior publication. However, we may already mention that 
cooling curves starting at low entropy can exhibit an initial increase 
in their luminosity due to deuterium burning. We subsequently no- 
ticed that the colder starts in fig. 8 of iMolliere & Mordasinil j2012h 
show a similar behaviour. 

To construct a model with a given mass M and internal 
entropy S, we integrate inwards from the photosphere and out- 
wards from the center, adjusting the central pressure, the cool- 
ing time t<j, and the luminosity L and the radius R at the top 
of the inner zone until the two integrations match at a pres- 
sure of 30 kbar. We use the Eddington approximation at the 
photosphere, setting T = T^g at Pphot = '^gl^K, and take the 
solar-metallicity (based on Lodders 2003) radiative opacity from 
iFreedman. Marlev. & Loddersl ( l2008l) . In the deep int erior, the con- 
tribution from the electron conductive opacity of ICassisi et alj 



We do make the standard assumption that the interior is fully convec- 
tive. even though stabilising compositional gradients have been suggested 
to possibly shut ofF large -scale convective movements ( iStevensonii 19791 : 
iLeconte & Chabrieiil2012l) . 
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( l2007h is also included. The equat ion of state for the hyd r ogen- 
helium mixture is that of Sa umon. Chabrier. & van HomI d 19951 
hereafter ISCvHI) with a helium mass fraction Y - 0.25. Since we 
focus on gas giants, we do not include a rocky core, which was 
found in a test grid to increase t he luminosities by at most a few 
percent at the lower masses, as in lSaumon et al. I (ll996h . 

The grid has a lower entropy limit of 5 - 7-7.5 ^B/baryorQ. 
The upper S limit in the grid is set by the requirement that R< 8 Rj , 
a value found to make numerical convergence straightforward. The 
upper limits are 5max - 12 near 1 Mj and 5 max - 14 for the higher 
masses. 

Fig. [T] shows interior profiles in the T-P plane for a range of 
entropies and masses. Schematically, since hydrostatic balance dic- 
tates ~ M-/R*, increasing the mass at a fixed entropy extends 
the centre to a higher pressure along the adiabat. This is exacer- 
bated at high entropies, where the planet substantially shrinks with 
increasing mass (while low-entropy objects have a roughly con- 
stant radius). At fixed mass, increasing the entropy mostly shifts 
the centre to higher T or to lower P. The first case obtains for 
low-entropy objects, which are essentially at zero temperature in 
the sense that k^T^ <c Ef, where is the central temperature 
and Ef the Fermi energy level at the centre. Increasing the en- 
tropy partially lifts the degeneracy since the degeneracy parameter 
y = k^T^IEf oc T^lg^-^^ is a monotonic function of S , and P^ re- 
mains constant because of the constant radius. At entropies higher 
than a turn-over value of 10 fes/baryon, the central temperature 
does not increase (and even decreases) with entropy. As for P^, it 
decreases because the radius increases. These behaviours also hold 
at higher masses not shown in the figure. 

As for the 1-Mj planet with 5=9 shown in Fig. [T] some 
models with entropy S - 8-9.5 show a second, 'detached' con- 
vective zone at lower pressures, which follows fro m a re-increase 
of the Rosseland mean opacity (see section 3.1 of ISurrows et alj 
fl997h . This second convection zone, which is at most at an en- 
tropy 0.2 fcs /baryon higher than the convective core, will not affect 
the evolution of the object since the radiative thermal timescale is 
much shorter than the cooling time ts , throughout the atmosphere. 
This holds in particular at the inner radiative-convective boundary 
(RGB); for example, the 1-Mj, 5 = 10 model has ts = 14 Myr and 
fthcnii = 0.2 Myr at its RGB. Finally, we note that the higher-entropy 
objects (5 > 11.5, with some dependence on mass) are convective 
from the centre all the way to the photosphere. 
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Figure 1. Temperature-pressure profiles of core-less planets with entropy 
5 = 7.5-12 A'B /baryon (bottom to top), in steps of 0.5 <:b /baryon, and 
M = 1,3, and 10 Mj (full, dashed, and dotted lines, respectively). The radii 
at5 = 7.5, 8, 9, 10, 11, and 12 /tB/baryon are 1.12, 1.18, 1.34, 1.66, 2.48, 
and 7.94 Rj (1 Mj), 1.11, 1.16, 1.29, 1.52, 1.96, and 2.85 Rj (3 Mj), and 
1.03, 1.07, 1.19, 1.40, 1.76, and 2.34 Rj (10 Mj), respectively. Squares indi- 
cate the centre and triangles mark where the enclosed mass is 99 per cent of 
the total mass for each model. Curves end at the photosphere in the Edding- 
ton approximation, and the thick part(s) of each curve show the radiative 
zone(s), while the thin part(s) of the curve are convective; the radiative- 
convective boundary relevant for the thermal evolution is highlighted by a 
ring. Thin grey lines at low P and T indicate contours of constant opac- 
ity iFreedman et al. 2008), with log*: = -3.5 to -1.0 (cm^ g"') in steps of 
0.5 dex, increasing with P. 

To try to understand these luminosity scalings, we firstly note 
that the radiative luminosity at the radiative-convective boundary, 
which is eq ual to the total luminosity, can always be written as 
jArras & B ildsten 2006) 

, AacT^dT] 167tGflc /Wr;'^,Vad,cb 

L = \4nr' = , (2) 

\ 3Kp dr /j.^^ 3 A-rcb/'rcb 

approximating the convection zone to contain the whole mass (cf. 
labels in Fig.[T]l and radius. Thus, one way of obtaining L{M,S) is 
to express the four quantities Pjcb, T^rcb, '^■jiAicb and ^Cpcb in terms of 
M and S . This will now be done for the three regimes regimes in 
turn, starting with low entropies. 



2.2 Luminosity as a function of mass and entropy 

To provide a model-independent way of thinking about an object's 
brightness and thus to facilitate comparison with other models, we 
show in Fig. [2] the luminosity of the planet as a function of its in- 
ternal entropy. We focus on objects without significant deuterium 
burning. Two main regimes are apparent: at lower entropies, the 
scaling with mass is roughly L ~ M, while at high entropies, the lu- 
minosity becomes almost insensitive to mass, i.e. L ~ M". Looking 
more closely, the high-luminosity regime is described by L oc M"^ 
(for M > 1 Mj), and the brief flattening of the luminosity slope 
with entropy at 5 ~ 8.2-9, depending on mass, marks a transition 
from L oc M' to L oc A/"^, respectively at fixed intermediate and 
low entropy. 

^ Entropy values in this work are always given in the usual units of mul- 
tiples of Boltzmann's constant per baryon (i.e. per mass of hydrogen 
atom). For comparison, 4.5-Gyr-old Jupiter currently has an entropy of 6.5. 



2.2.1 Low-entropy regime: L oc M^'' fi^wiS ) 

To begin, consider the entropy dependence of luminosity at a fixed 
mass when 5 < 9-8 (at 1-10 Mj respectively). Fig.[T]reveals that 
the opacity at the radiative-convective boundary A-^b remains con- 
stant at a given mass, which provides a first relation. Secondly, over 
the temperature and pressure range of interest in this regime, the 
hydrogen and helium remain respectively molecular and neutral. 
This implies that Vadicb is constant and that the entropy has a sim- 
ple functional form in the idea l-gas approxi mation, given by the 
Sackur-Tetrode expression (e.g. lGallenll985h . For an H2-He i mix- 
ture with y = 7 , this is 

4 

S = 9.6 + ^ ln(r/1600 K) - \n(P/3 bar) (3) 
32 16 

(hence Vad = 0.31), where In 10 = 2.3, and the reference T and P 
values were chosen for Section l2.2.2l Finally, r^cb is approximately 
fit by r,cb oc 10" "* . Gombining the four relations (constant V^drcb 
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9 10 
Entropy (ka/boryon) 

Figure 2. Luminosity as a function of entropy for our planet models. The 
masses shown are (bottom to top) M = 1, 3, 5, and 10 Mj. Here as through- 
out this paper, we use Lsun = 3.86 X 10^^ erg s"' . 



and the three non-trivial ones) with equation (|2j yields L oc 10'"''* . 
This is quite close to a power-law fit of Fig. |2l which gives L oc 
10' 

Next, consider the position of the RGB at fixed low entropy for 
different masses. Fixed entropy immediately implies T oc P^'iJ, and 
the constant V^d = 0.31 gives the second relation. The profiles of 
Fig.[T]indicate that in this low-entropy regime, the opacity depends 
approximately only on the pressure, with k oc P" '. Finally, we find 
that Prcb oc Af-^, thus providing the fourth relation. Combining 
these and using equation l|2j gives L oc M"^, in good agreement 
with a direct fit, which gives L oc M" ' . 

Thus, a constant adiabatic gradient, an RCB opacity depen- 
dent only on the mass, and a few power laws suffice to show that in 
the low-entropy regime 



L,™ = 10-"Lo(^ 



10 



1.3(5-7.5) 



(4) 



which defines fiow(S) up to a constant. For Jupiter's adiabat with 
T = 165-170 K at 1 bar jSaumon & Guillolll2004) and thus an 
entropy of 6.71-6.75 /rB/baryori, this predicts =^ 1.9 x 10"' Lq, in 
good agreement with the current value Lj = 8.7 x 10"'" Lq. 



2.2.2 Intermediate-entropy regime: L oc Mf^iS) 

We now look at intermediate entropies, which are in the range 9- 
10.2 (at 1 Mj) or 8.2-9 .5 yts/baryon (at 10 Mj). As explained by 
lArras & BildstenI ( l2006h . equation l|2j is of the form L = Mf(S ), 
where f(S ) is a function of the entropy in the convection zone, 
only if the quantity T'^V-^ci/kP at the RCB is a function of a unique 
variable, S . The intermediate-S behaviour can then be understood 
by noticing in Fig. [T] that objects at those entropies have a r,cb 
independent of internal entropy and an extended atmosphere (in- 
terrupted or not by a second convection zone). For these planets, 
the photosphere is sufficiently far from the RCB (Pjcb -SJ lOPphot, 
Ticb'* ^ 5 T'eff'*) that t he atmosphere merg es onto the radiative-zero 



profile, which solves jCox & Giulilll968 ) 
dT 3 L K 
dP 



(5) 



16acG M T^' 

where the boundary condition is by definition of no consequence. 
The solution is a T(P) relation, which yields the radiative gradi- 
ent Vrad = dlnr/dln/* along it. The important point is that both 



the atmosphere profile and its gradient depend only on the quantity 
L/M. Now, choosing an internal entropy fixes the adiabat, i.e. sets 
a second T(P) relation. We require that at the intersection of the 
atmosphere and the adiabat Viad be equal to Vaj, which is the slope 
of the chosen adiabat. This thus pins down Vi-^d in the atmosphere 
and hence L/M (and the atmosphere profile itself). Therefore, there 
is a unique L/M associated to an 5 , which means that L must be of 
thefom0L/M = /(5). 

Since we use the Schwarzschild criterion V;,d < Vmd to deter- 
mine convective instability, where the radiative gradient is 



V,-ad 



3LP 



\6nacGMT* 



(6) 



near the surface, and Vad is relatively constant, a slow inward in- 
crease of K will ensure a deep RCB. Consequently, the L oc Mf(S) 
scaling will hold for those M and 5 such that, starting at the pho- 
tosphere with Pphot - f{M) and T^ff ^ f(Pphot,S), the opacity in- 
creases only slowly along the adiabats. (For a power-law opacity 
K = K(,P"T^"^\ this means n/(n + s) ^ V^d-) This is indeed the case 
for intermediate-S models in Fig. [T] where the T(P) profiles are 
nearly aligned with contour lines of constant opacity. For this ar- 
gu ment to hold, the radius mu st be rather independent of the mass 
(cf. lZapolskv & Salpeteill96^ and of the entropy. Also, one needs 
to approximate the entropy at the photosphere to be that of the in- 
terior, which is reasonable: even in the extended atmospheres, the 
entropy increases by at most ^ 0.5 A;B/baryon. 

The functional form of /(5) can be obtained by fixing Tpcb 
and using equation ^ for the entropy of an ideal gas. The refer- 
ence temperature of 1600 K was chosen based on the radiative- 
convective boundaries of Fig. [T] For convenience, 5(3 bar, 1600 K) 
was computed using the interpolated SCvH tables. These include 
the contribution from the ideal entropy of mixing a remarkably 
constant 5 mix = 0.18 A:B/baryon in a large region away from 
H2 dissociation. Combining this S (P, T) with equation Q, fixing 
T = 1600 K, and taking k(P, 1600 K) ^ 0.0104(/'/3 bar)"' cm^ g"' 
from the interpolated table yields 



logi 



M/Mj 



-5.05 + 1.51(5 -9.6), 



(7) 



ergs-'Mj-' lO'-SX^-w). 



i.e. L = Mfy-AS) with /„(5) = 3.4 x lO^** 
(Thus, ifKoc P" and d5/d lnP = -b, log,o MS ) oc (« + 1)/^? In 10.) 
The subscript 'rz' highlights that the solution applies when the ra- 
diative zero solution is reached. This excellently fits the L/M-S 
relationship found in our grid of models at intermediate entropies. 

2.2.3 High-entropy regime: L oc M" '/hii,h(5) 

For 5 > 10-9 (at 3-10 Mj, respectively), the luminosity becomes 
almost independent of mass at a given entropy. This indicates that 
the radiative solution does not hold anymore, and indeed Fig. [T] 
shows that planets with high entropy have atmospheres extending 
only over a small pressure range, with the more massive objects 
fully convective from the centre to the photosphere. The shortness 
of the atmosphere is due to the opacity's rapid increase inward, 
as constant-/c contours are almost perpendicular to T(P) profiles in 



This regime was obtained bv lArras & BildstenI b006l) when looking at 
irradiated planets. This can be roughly thought to fix Tr^b to the irradiation 
temperature, so that T^V.^^/kP at the RCB is automatically a function of 
only one thermodynamic variable, for instance 5, for all e ntropies. 

We used the con'ected vers ion of the equations in ISCvHt see the appendix 
of lSaumon & MarlevI i2003) . 
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that region. As in Section 12.2.1 1 we look at the behaviour of T, P, 
K, and V;,(] at the RCB as a function of M and S . 

Fig.[T]shows that at a fixed mass, r^b is almost indepedent of 
the (high) entropy, with the actual scaling closer to T^cb* 10""'. 
Also, ATrcb is again mostly independent of 5 , as for the low entropies. 
At high entropies, Vi,ci,.;.|, drops continuously with increasing S (de- 
creasing Prcb), with a very rough Vaj oc lO"" ""*. For the fourth 
relation, we can fit Prcb lO""^'' . Combining all this with equation 
(|2) then gives L oc lO"^'' , which is quite close to a fit L oc lO"'^'^ . 

At fixed high entropy, Vad is somewhat constant at 0.25-0.13 
for S = 10.5-12. Also, above ^ 1500 K and at P < 1 bar, the 
opacity is approximately independent of pressure, scaling only with 
temperature as /c ~ T*. Finally, as in the low-entropy regime, Pjcb oc 
M"-'; combining the four relations, we should have L oc M" with 
a ^ 0.5-0.4. This is not far from the direct fit L oc M"^. 

Therefore, the luminosity in the high-entropy regime is 



10' 



,0.6(5-10.5) 



(8) 



which defines /high (5) up to a multiplicative constant. Fig. |2] indi- 
cates that the coefficients are slightly different at the lowest masses 
(M < 1.3 Mj as it turns out), with the luminosity more sensitive on 
the mass, but this was ignored when obtaining equation ([8}. 

Before summarising, let us briefly digress about the Pjcb oc 
M"-^ scaling seen both at fixed low entropy and at fixed high en- 
tropy. In both cases, the model grid shows that, as a reasonable 
approximation, P^cb oc Pphot- At low entropies, R is constant, such 
that Prcb fphot = '^gl^K °^ M/k. Then, the found k ~ P (see 
Section l2.2. It immediately implies Pj^b °^ M"^. At high entropies, 
planetary radii are significantly larger and vary, such that the radius 
dependence of the photospheric pressure should not be neglected; 
thus Pphoi oc M/R^K. Fitting the R(M, fixed high 5) relation in our 
models, we find R ~ M"' ' or R ~ M"'^, depending on the entropy. 
Then, ^ 0.2 and k ~ T* yield k ~ P"^ at fixed 5, and thus, 
isolating, P,cb ~ or P,cb ~ M" '. This is a slightly stronger de- 
pendence on mass than what is found in the grid, but the argument 
shows how the rough scaling can be derived. 



2.2.4 Summary of luminosity scalings 

In summary, we found from fitting the L(M,S) relation across 
our planet models that L oc M" 10'"', with a = (0.7, 1,0.3) and 
A = (1.3,1.5,0.6) at low (S < 9-8, for 1-10 Mj respectively), 
intermediate (5 < 10.2-9.5) or high entropy, respectively. These 
scalings can be fit directly or more or less accurately derived by 
considering the behaviour of the different factors in equation (|2)- 

Our approximate understanding of the different regimes is 
thus the following. For the conditions found in the atmospheres of 
intermediate-entropy planets, contours of constant opacity are al- 
most parallel to adiabats. Since Vrad oc /<• by equation |6] and, for a 
constant V^d, only Vjad determines when the atmosphere becomes 
convective, a slow inward increase of the opacity causes the radia- 
tive zone to extend over a large pressure range. This in turn means 
that the atmosphere can reach the radiative-zero solution, which 
we have shown necessarily implies L = Mf(S). For the condi- 
tions found in the atmospheres of planets with high and low en- 
tropy, however, opacity increases relatively quickly along an adi- 
abat. Since the T(P) slope in the atmosphere is not too different 
from that of the convective zone's adiabat, this means that opac- 
ity increases rapidly in the atmosphere, which therefore cannot 
join onto the radiative-zero solution before becoming convective. 
The diff'erent L{M,S) scalings then reflect in part the approximate 



temperature- and pressure-independence of the opacity at low and 
high 5 , respectively. It is interesting to note that the transition from 
low to intermediate entropies is accompanied by a 'second-order' 
(i.e. relatively small) change in a and A, while the the physical ex- 
planation changes to 'zeroth order'. 



2.3 Cooling 

With L oc M" 10-'^ equation Q implies d5/df = -M"-"^ f{S)IT . 
Since hydrostatic balance yields P^ ~ GM^/R* and the convective 
core is adiabatic, we expect T ~ oc {M^ jR^ f"^ . (Across the 
whole grid, T IT^ = 0.55-0.63.) Therefore, at constant entropy and 
ignoring the radius dependence, -d5/d/ oc M"^^'"'"'^". This means 
that the entropy of a cooling planet should be approximately a func- 
tion of ?/M^'-"i when a ^ 1 (at low and intermediate entropies), 
and that more massive planets cool more slowly (since a is always 
smaller than 2Vad 1 - 5/3). Thus, at fixed entropy, L oc M", 
with a = 0.3-1, but at a fixed time the luminosity has a steeper de- 
pendence on M, closer to L oc M", because more massive planets 
move to lower entropy more slowly. A more detailed analytic un- 
derstanding of the coolin g curves for irradiated planets along these 
lines was developed bv lArras & BildstenI ( |2006|) . 



2.4 Comparison with previous work 

We now compare our cooling curves to classical hot starts. Fig. |3] 
shows cooling curves^ for large initial_entropy compared to the 
hot-start models of iMarlev et al.l ilOOlj). the COND03 models 
tearaffe et all l2003b and those of iBurrows et all ( Il997h . all of 
which use non-grey atmospheres with detailed opacities. The 
agreement is excellent, with our luminosities within the first 3 Gyr 
approximately within -30 an d 2 per cent and - 30 and 60 per cent 
above those of lBaraffe et alj and lBurrows et al.l (Of interest for the 
example of Section O our radii along the cooling sequence are at 
most approximately 0.08 an d 0.04 Rj (8 and 4 per cen t) greater, 
respectively, as mentioned in iHuang & Cumminal ( l2012l) . This dif- 
ference is comparable to the effect of neglecting hea vy elements in 
the eq uation of state or not including a solid core jSaumon et alj 
1 19961) and not significant for our puiposes.) Our deuterium-burning 
phase at 20 Mj ends slightly earlier than in iBurrows et al.l d 19971) 
but this might be due to our simplified treatment of the screening 
factor 

Fig.|4]shows cooling curves for lower initial entropi es tha n in 
Fig. [3] The cooling curves show the behaviour found by "mot" and 
Spieg el & Burrow si ( I2OI2I) in which the luminosity initially varies 
very slowly, with the cooling time at the initial entropy much larger 
than the age of the planet. Eventually, the cooling curve joins the 
hot start cooling curve once the cooling time becomes comparable 
to the age. 

Comparing to the cold starts in fig. 4 of lM07l our models 
are a factor of - 3.5-3.9 lower in luminosity for the same initial 
entropies, as given in their fig. 2. Increasing our initial entropies 
by 0.38 fee/batyon brings our cooling curves into agreement with 
theirs when the planet has not yet starting cooling. (For this com- 
parison, we do not c orrect the time offset for the higher masses 
(see figs. 2 and 4 of iMOTh . which are already on hot-start cool- 
ing curves at the earliest times shown.) Ispiegel & Burrows] (l2012h 
computed the evolution of gas giants starting with a wide range 
of initial entropies. To compute the bolometric luminosity of their 
models, we take the published spectra and integrate the flux in the 
wavelength range given, 0.8-15 vim. The bottom panel of Fig. |4] 
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Figure 3. Cooling curves for (bottom t o top) M = 1 , 2, 3, 5 , 8, 10, 15, and 
20 Mj (solid cui'ves) compared with the lMarlev et al.lfcoOTh hot starts (for 
M = 1, 2, and 10 Mj o nly; dash-dotted curves), the COND03 tracks (dots) 
and the lBurrows et alj (il997) models (dashed curves). The COND03 data 
is shown as dots because of insufficient sampling of the cooHng tracks at 
higher masses. 



shows the comparison to all four models, with or without clouds 
and at solar or three times solar metallicity. Increasing our en- 
tropy by 0.40 ^a/baryon - e.g. comparing the ISpiegel & Burrowsl 
( I2OI2I) model with 5 = 9.0 to ours with 5 - 9A - yields very 
good agreement, with our luminosities overlying their curves or 
within the spread due to the different atmospheres. The apparent 
disagreement at late times comes from the increasing fraction of 
the flux in the Ray leigh- Jeans tail of the spectrum beyond 15 |J.m. 
For comparison, the implied required bolometric correction is equal 
to 10 to 50 per cent of the flux in 0.8-15 (im for a blackbody with 
TeB - 700-300 K. From t his and the hot-start T^ff tracks shown in 
ISpiegel & BurrowH (1201 2h . we estimate that integrating the spec- 
trum should give a reasonable estimate (to ca. 30 per cent) of the 
bolometric luminosity only up to ^ 50, 200, and 1000 Myr for ob- 
jects with M = 1, 3, and 10 Mj respectively. This is indeed seen in 
Fig.! 

In contrast with this is our nearly perfect agreeme nt with the 
MESA stellar evolution code ( Paxton et alj|20lil l2013h . We com- 
pared our L{S) relation to ones obtained from mesa at different 
masses (as in Fig. O and found that they are very nearly the same, 
with an entropy offset AS < 0.05. We also produced grids with 
other opacities, using the mesa tables with the default Y = 0.28 
and Z = 0.02, and the lFreedman et al.l ( |2008|) tables with [M/H] = 
±0.3 dex, and found that this changed the luminosity by at most 
=: 10 per cent. Similar differences were found to result from a 
changed helium mass fraction in the bulk of the planet. Therefore, 
it is not clear what might explain the small disagreement between 
the entropy values of these various au thors. O ne possibility is a dif- 
ference in the implementation of the ISCvHI equation of state, but 
this remains to be investigated. 

The upshot of these comparisons to classical, non-grey- 
atmosphere hot and cold starts is that cooling tracks computed with 
the simple and numerically swift cooling approach described above 
can reproduce models which explicitly calculate the time depen- 
dence of the luminosity. When comparing models from different 
groups, one should keep in mind that there is a systematic uncer- 
tainty in the entropy values of A5syst - 0.4. This is however small 
compared to the entropy difference between hot and cold starts. 
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Figure 4. Top pan el: Cooling curves f or M = 1, 2, 4, 10 Mj and initial en- 
tropies equal to the lMarlev et al.l fcoOVl MOV) 'tuning-fork diagram' values 
5i = 9.23, 9.00, 8.60, 8.23 (see thei r fig. 2) increased by 0.38. The dotted 
curves show the lMarlev et al.l bOOVh cooling tracks. Bottom panel: Cooling 
curves for M = 1,3, 10 Mj (bottom to top) with initial entropies Si = 9.0 
and 5 i = 1 1 incre a sed b y 0.40 (solid lines), compared to the models of 
ISpiegel & BuiTowsl | |20121 SB12) with S, = 9.0 and S, = II, to all four 
atmosphere types at each mass (dotted; hybrid-cloud or cloud-free, solar 
or thre e times solar metallicity). The truncation of the ISpiegel & Burrowsl 
l' 2012l) spectra causes an underestimate of the bolometric luminosity at late 
times (see text). 



3 GENERAL CONSTRAINTS FROM LUMINOSITY 
MEASUREMENTS 

Masses of directly-detected exo planets are usually inferred by fit- 
ting h ot- start cooling curves (e.g. lBurrows et al .ll997l;lBaraffe et alj 
l2003h to the measured luminosity of the planet, using the stellar age 
as the cooling time. Since the hot-start luminosity at a given age is a 
function only of planet mass, the measured luminosity determines 
the planet mass. Moreover, a planet's luminosity, at a given time 
and for a given mass, can never exceed that of the hot starts, since 
a larger initial entropy would have merely cooled onto the hot-start 
cooling track at an earlier age. 

However, we have seen above that the luminosity at a given 
mass can be lowered by considering a sufficiently smaller initial 
entropy, which might be the outcome of more realistic formation 
scenarios ( iMarlev et aT]|2007l; ISpiegel & BurrowsllioTih . With the 
fact that luminosity increases with planet mass, this simple state- 
ment has important consequences for the interpretation of direct- 
detection measurements, namely that there is not a unique mass 
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which has a given himinosity at a given age. Cold-start solutions 
correspond to planets not having forgotten their initial conditions, 
specifically their initial entropy S \, and every different initial en- 
tropy is associated to a different mass. In other words, a point 
in {t, L) space - a single bright ness measurement - is mapped to 
a curve in (M, 5i) space. Since iMarlev et al] (12007^ (but see also 
lBaraffeetal.ll2002l: l Fortnev et aLlllOOSh . it is generally recognised 
that direct detections should not be inter preted to yield a uniqu e 
mass solution, but, with the exception of iBonnefov et alj j2013h . 
who used infrared photometry, this is the first time that this degen- 
eracy is calculated explicitly. 



3.1 Shape of the MS-, constraints 

The top panel of Fig. |5] shows the allowed masses and initial en- 
tropies for different values of luminosity Li,oi/Lq = 10"^, 10"* 
and 10"^, at ages of 10, 30, and 100 Myr. Below the deuterium- 
burning mass, constant-luminosity curves in the M-S\ plane have 
two regimes. At high initial entropies, the derived mass is the 'hot- 
start mass' independent of Si since all Si greater than a certain 
value have cooling times shorter than the age of the system. There, 
uncertainty in the stellar age translates directly into uncertainty in 
the planet mass: since L oc 1// and L ~ M', the mass uncertainty 
is AM/M ^At/t. At lower entropies, the luminosity measure- 
ment occurs during the early, almost constant-luminosity evolution 
phase. Given a luminosity match in this region, one can obtain an- 
other by assuming a lower (higher) initial entropy and compensat- 
ing by increasing (decreasing) the mass. As seen in Section |2] L 
is a very sensitive function of S at low and intermediate entropies, 
so that a small decrease in initial entropy must be compensated by 
a large increase in mass to yield the same luminosity at a given 
time; this yields the approximately flat portion of the curves. As 
long as the cooling time for a range of masses and entropies re- 
mains shorter than the age, the entropy constraints do not signifi- 
cantly depend on the age. The uncertainty in the initial entropy is 
ASi - 1//1 AlogigLboh where l/A - 0.7 or 1.7 at low or high en- 
tropy (see Section [2.2.4) . Finally, at higher masses, the possibility 
of an initial increase in the luminosity makes L(t) non-monotonic 
and hence results in more complex geometries in the M-Si plane. 
This will be discussed in a forthcoming publication. 

The circles in the top panel of Fig.|5]show the initial entropies 
for cold- and hot-start models from lM07L the 'tuning-fork diagram 
of exoplanets'. The entropies were increased by 0.38 ^s/baryon as 
in the top panel of Fig. |4]to match our models. Since the cooling 
time increases with mass (see Section |2j3}, heavier planets of the 
hot-start (upper) branch, i.e. with arbitrarily high initial luminos- 
ity, have cooled less and are therefore at higher S . For their part, 
the cold-start entropies, which are still the post-formation ones, lie 
close to a curve of constant luminosity L ^ 2 x 10"* Lq. This re- 
flects the fact that the post-formation luminosities in lMOTl as seen 
in their fig. 3, have similar values for all masses. The cause of this 
'coincidence' is presently not clear; it might be a physical process 
or an artifact of the procedure used to form planets of different 
masses. Putting uncertainties in the precise values aside, the two 
prongs of the tuning fork in Fig. |5] give an approximate bracket 
within which just-formed planets might be found. 

Mass information for a directly-detected planet can put useful 
constraints on its initial entropy and also potentially on its age and 
luminosity simultaneously. For instance, dynamical-stability anal- 
yses and radial-velocity observations (see Section|5} typically pro- 
vide upper bounds on the masses. Since Si decreases monotonically 
with M at a fixed luminosity, this translates into a lower limit on the 
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Figure 5. Top panel: Allowed values of planet mass M and initial en- 
tropy Si corresponding to Lboi = 10"', 10"*, and 10"' (bottom to top 
sets of curves) and ages t = 10. 30, and 100 Myr (left to right or dot- 
ted, full and dashed lines). The circles show the results from lMarlev et alj 
(2007), increased by 0.38 /cfl/baryon (cf. Fig.|4), for cold starts (filled cir- 
cles) and hot starts (open circles), where the entropy is 1 Myr after the 
onset of cooling. Red and blue squares indicate the illustrative M and Si 
used for the lightcurves in the bottom panel, with errorbars correspond- 
ing to the hatched regions there. Bottom panel: Two examples of cooling 
curves that have L = 10"* Lq at 30 Myr The dot-dashed (red) line has 
M = 1.85 Mj and Si = 11.9 (hot start), while the full (blue) one has 
M = 12 Mj and Si = 8.33 (cold start). Hashed regions indicate a spread 
of 0.5 Mj (0. 1 /baryon) for the hot- (cold-)start curve. The respective rel- 
ative uncertainties are AM/M =: jAt/t and ASi - O.VAlogjQ Lboi for the 
cold- and hot-stait cases, respectively (see text). 



initial entropy. This has the potential of excluding the coldest-start 
formation scenarios. Conversely, a lower limit on the planet's mass 
implies an upper bound on Si. If it is greater than the hot-start mass, 
this lower limit on the mass of a planet would be very powerful, due 
to the verticalness of the hot-star branch. Combined with the flat- 
ness of the 'cold branch' of the M(Si) curve, this could easily re- 
strict the initial entropy to a dramatically small ASi - 0.5 ~ ASsys(. 
Also, the top panel of Fig. [5] shows that not all age and luminosity 
combinations are consistent with a given mass upper limit. Given 
the often important uncertainties in the age and the bolometric lu- 
minosity, this may represent valuable input. 
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Figure 6. Entropy of hot staits as a function of time for planet masses of 1, 
3, and 10 Mj (bottom to top). At a given age, the curve indicates the value of 
initial entropy above which the 'hot-start mass' applies. For a planet mass 
larger than the hot-start value, the initial entropy must be lower than the 
hot-start entropy at the cun'ent age. 



3.2 Solutions on the hot- vs. cold-start branch 

The bottom panel of Fig. [5] shows lightcui-ves illustrating the two 
regimes of the M(Si) curves discussed above. The hot-start mass is 
1.85 Mj, whereas the selected cold-start case (5i = 8.33) has M = 
12 Mj - i.e. a six times larger mass - and both reach log Lboi/LQ = 
-6 at f = 30 Myr, with the cold-start values essentially independent 
of age. The hatched region around the 1.85-M] curve comes from 
hot-start solutions between 1.35 and 2.35 Mj and is within a factor 
of two of the target luminosity, showing the moderate sensitivity of 
the cooling curves to the mass. However, in the cold-start phase, 
a variation by a factor of two can also be obtained by varying at 
a fixed mass the initial entropy from 8.23 to 8.43 /cB/baryori. This 
great sensitivity implies that combining a luminosity measurement 
with information on the mass would yield, if some of the hot-start 
masses can be excluded, tight constraints on the initial entropy. 

3.3 Deflnition of 'hot-start mass' 

By showing the entropy of hot starts as a function of time, Fig.|6] 
provides another way of looking at Fig. |5] Given a mass obtained 
from hot-start cooling tracks, the entropy value at that time read 
off from the curves indicates what 'hot' is, i.e. provides a lower 
bound on the initial entropy if the hot-start mass is the true mass. 
If however the planet is more massive, this entropy value is instead 
an upper bound on the post-formation entropy. 

Since it has no bearing upon the detections considered be- 
low, we defer a similar discussion concerning deuterium-burning 
objects to a later work. As noted before, the possibility of an initial 
increase in the luminosity with time can complicate the form of the 
M(Si) constraints. 



gravity by fitting its photometry and spectra with atmosphere mod- 
els. Integrating the best-fitting model spectrum gives the luminos- 
ity, and this can be combined with T^ff and log g to yield the radius 
and the mass. This proced ure w as described and carefully applied 
by_ Mohantv et alj ( 12004) and iMohantv, Javawardhana. & Basril 
12004'). However, one can go further: considering models coupling 
detailed atmospheres with interiors at an arbitrary entropy, the mass 
and radius translate into a mass and current entropy. Then, using 
cooling tracks beginning with a range of initial conditions, the ini- 
tial state of the object can be deduced given the age. Thus, in con- 
trast to the case when only luminosity is used, M and 5 i can both 
be determined without any degeneracy between the two. 

In practice, however, there seems to be too much uncertainty 
in atmosphere models for this method to be currently reliable, as 
the work of Mohantv et al. shows. Their sample comprised a dozen 
young (~ 5 Myr) objects with T.s ^ 2600-2900 kH with high- 
resolution optical spectra. The combined presence of a gravity- 
sensitive Na I doublet and effective-temperature-sensitive TiO band 
near 0.8 fi.m allowed a relatively precise determination of log g and 
reff for most objects , with uncertainties of 0.25 dex and 150 K 
jMohantv et"ai]l2004h . However, there were significant offsets in 
the gravity (- 0.5 dex) of the two c o ldest object s with respect 
to model predictions of Baraffe et alj ( 1 19981) and IChabrier et al.l 
( l2000h . lMohantv et alj came to the conclusion that the mod- 

els' treatment of deuterium burning, convectior0 or accretion - i.e. 
the assumed initial conditions - are most likely responsible for 
thi s disagreemen t at low er T^f;. Moreover, the more recent work 
of lBarmanetallfeOllbl) (see Section 1 . 1 1 below) indicates that 
'unexpected' cloud thickness and non-equilibrium chemistry may 
compromise a straightforward intepretations of spectra in terms of 
gravity and temperature for young, low-mass objects. Nevertheless, 
with the hope that future observations will allow a reliable calibra- 
tion of atmosphere models, we illustrate with an example how M 
and can be determined for an object from its measured T^ff and 

Fig. |7] shows the constraints on the current mass and entropy 
of a planet with log g = 4.00 ± 0.25 (cm s^^) and Jeff = 900 ± 50 K, 
where the uncertainties in the gravity and effective t emper ature 
are the possible accuracy reported by iMohantv et all (|2004 and 
thus correspond to an optimistic scenario. The constraints are ob- 
tained by simultaneously solving for M, S , R, and L given T^ff 
and log^, with the L(M,S ) and R(M,S ) relations given by our 
model grid. The required mass and current entropy are M = 6.5 Mj 
and 5 = 9.2, with the l-cr ellipse between 4.0-10.1 Mj and 9.0- 
9.3 ^g/batyon. The large uncertainty in the mass (the width of the 
ellipse) is dominated by the uncertainty in the gravity, since radius 
is roughly constant at these entropies. We note that non-Gau6ian er- 
rorbars are trivial to propagate when determining mass and entropy 
in this way since it is only a matter of mapping each (log g, T^g) 
pair to an (M, S ) point. 

Since this MS determination concerns only the cutTent state, 
it is independent of the cooling sequence, in particular of the 'hot 
vs. cold start' issue. Nevertheless, with this approach, it is immedi- 



4 GENERAL CONSTRAINTS FROM GRAVITY AND 
EFFECTIVE-TEMPERATURE MEASUREMENTS 

Before applying the analysis described in the previous section to 
observed systems, it is worth discussing a second way by which 
constraints on the mass and initial entropy can be obtained. The 
idea is to firstly derive an object's effective temperature and surface 



^ The spectral types are =: M-M7.5, but lMohantv. Javawardhana, & Basril 
stress that the correspondence between spectral type and effective 
temperature of young objects has not yet been empirically established and 
thus that calibration work (in continuation of theirs) remains to be done. 
* Qualitatively, their finding that theoretical tracks predict too quick 
cooling at low masses might be evidence for the argument of 
iLeconte & Chabriej Jioll) that convection in the interior of these objects 
could be less efficient than usually assumed. 
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Figure 7. Constraints on the mass and entropy of an example planet from 
its gravity and effective temperature, with log g = 4.00 ± 0.25 (cm s"~) and 
= 900 ± 50 K (and thus optimistic errorbars; see text). The mass and 
current entropy corresponding to log g = 4.00 and T^g = 900 K are shown 
by a filled circle. The narrow ellipse (solid line) is the l-o" confidence region 
in the mass and current entropy, whereas the deformed ellipse (dotted line) 
shows the constraints on the initial entropy assuming an age of 30 Myr The 
open circle corresponds to the 'measured' central value. The three solid 
grey lines indicate the entropy of hot starts at 20, 30, and 50 Myr (top to 
bottom). Higher masses can be seen to cool more slowly (see Section[23j- 



ately apparent what constraints the age imposes on the initial con- 
ditions. Indeed, not all (M, 5) are consistent with an age since no 
planet of a given mass can be at a higher entropy than the hot-start 
model at that time (see Fig.[6] The entropy of hot starts is shown in 
Fig. Rafter 10, 30, and 50 Myr. The 30-Myr age excludes objects 
with M < 6 Mj, which is the hot-start mass of this example. Con- 
sidering only hot-start evolution sequences would have been equiv- 
alent to requiring the solution to be on one of the hot-start (grey) 
curves. This is however a restriction which currently could not be 
justified, given our ignorance about the outcome of the formation 
processes. 

The dotted lines in Fig. |7] indicate the derived constraints on 
the initial entropy, assuming an age of 30 Myr These constraints 
are similar to ones based on luminosity (see Fig. |5} but are some- 
what tighter since an upper mass limit is provided by the measure- 
ment of logg. Even within a set of models, i.e. putting aside pos- 
sible systematic issues with the atmospheres of young objects, it 
is however often the case that the surface gravity is rather ill de- 
termined (as for the objects discussed below in Section |5}- In this 
case, provided a sufficiently large portion of the spectrum is cov- 
ered, we expect the approach based on the bolometric luminosity 
presented in this work to be more robust than the derivation of con- 
straints from effectively only the surface temperature. Indeed, the 
former avoids compounding uncertainties in T^s with those in the 
radius in an evolutionary sequence, which can be further affected 
by the presence of a core (of unknown mass). With both sustained 
modelling efforts and the detailed characterisation of an increasing 
number of detections, one may hope that reliable atmosphere mod- 
els for young objects will become available in a near future, allow- 
ing accurate determinations of the mass, radius, and initial entropy 
of directly-detected exoplanets. 



5 COMPARISON WITH OBSERVED OBJECTS 
5.1 Directly-detected objects 

iNeuhauser & Schmid ^ ( l2012h recently compiled and homoge- 
nously analysed photometric and spectral data for directly-imaged 
objects and candidates, selecting only those for which a mass be- 
low 25 M] is possibl^!]. They report for each of the 28 objects a 
luminosity and effective temperature, which they either take from 
the original papers or calculate, usually from bolometric correc- 
tions when a spectral type or colour in dex is available or bright ness 
difference with the primary when not. INeuhauser & Schmidt! then 
use a number of hot-start cooling models to derive (hot-start) mass 
values along with errorbars, while recognising that hot starts suffer 
from uncertainties at early ages. 

In this Section, we determine joint constraints on the masses 
and initial entropies of directly-detected objects, focussing on the 
ones for which (tentative) additional mass information is available. 
This is the n a generalisation, for the ob jects we consider, of the 
analysis of Neuhiiuser & S chtnidd ( 12012 '). A proper statistical anal- 
ysis of the set of M(S \) curves would be challenging at this point 
due to the inhomogeneity of the observational campaign designs. 
However, upcoming surveys should produce sets of observations 
with well-understood and homogenous biases, convenient for a sta- 
tistical treatment. 

Before turning to specific obj ects, we display in Fig . [8] the 
dat<E| collected and computed by INeuhauser & Schmidt! ( 12012!) 
along with hot- and cold-start cooling tracks for different masses . 
This is an update of th e analo gous figures of iMarlev et al.l ( 12007! , 
ImOtI) and lJanson et alj ( 1201 ih , which had only a handful of data 
points. This shows that there are already many data points which 

- at least based solely on their luminosity - could be explained 
by cold, warm, or hot starts, highlighting the importance of being 
open-minded about the initial entropy when inteipreting these ob- 
servations. 

There are two features of the data distribution in Fig. |8]which 
immediately stand out. The first is that the faintest young objects 
are brighter than the faintest oldest objects, i.e. that the minimum 
detected luminosity decreases with the age of the object. More- 
over, this minimum, with the exception of data points (7) and (28), 
approximately follows the cooling track of hot-start objects of 

- 7 Mj. The second feature is the absence of detections between 
the hot-start 10- and I5-Mj cooling curves, roughly between 20 and 
100 Myr More accurately, there is, in a given age bin in that range, 
a lower density of data points with luminosities around 10 Lq 
than at higher or lower luminosities. 

A proper assesment of the statistical significance of this 'gap' 
in the data points would require taking both the smallness of the 
number of detections around 40 Myr and the biases and non- 
detections of the various surveys into account. However, it would 
not be surprising if the underdensity in the luminosity function 
dN/dL of the data points, at a fixed age, proved to be real, since 
there is also a suggestive underdensity in the cooling curves. In- 
deed, the onset of deuterium burning near 13.6 Mj slows down the 
cooling, which breaks the hot-start scaling L oc r' (see Section [23] l 



' This value was chosen by ISchneider et al.l J201 ll) as an approximation 
to the 'brown-dwarf desert', which is a gap in the mass spectrum be- 
tween =^ 25-90 Mj (Marcy & Butler 200^; lOrether & Lineweave j liool 

Luhman et al. 2007; Dieterich et al. 2012). 

^ De spite doubts about its nature or existence (e.g. Jan son et alj |2012| ; 
iCurrie et al.ll20123) . we leave the re ported values for Fomalhaut b in this 
plot and also show the revised age of lMamaieS i2012l) . 
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Figure 8. Directly-observed objects with a hot-start mass < 25 Mj 
jNeuhauser & Schmidll r2012h compared with coohng curves for M = 
1, 2, 3, 5, 10, 15, 20, and 25 Mj (bottom to top), for 5i = 9.5 (dashed) and 
hot starts (dash-dotted); the 25-Mj curve is from lBurrows et al and 
has no cold-start equivalent. Since their uncertainties ai'c large, the central 
age values are taken as the geometric mean of the upper and lower bounds 
reported (only RXJ1609 B/b, GSC 8047 B/b, and HIP 78530 B/b have a 
reported best age, for all three of 5 Myr). The errorbars for the H R 8799 
planets go up to 160 Myr and do not include the lMova et aljfioiObh astero- 
seismology measurement of 1 . 1-1 .6 Gyr since it is not used in our analysis 
(see Section P^2.U . However, wider paler er rorbars with a ring symbol d o 
use the upper value of 1 . 1 Gyr (as reported bv'Neuh auser & Schmidll2012h . 
Note that data points (14) and (15) overlap, differing only in their verti- 
cal errorbai's. The upper limit for Fomalhaut b is included, also with the 
revised, greater age from Mamajek (2012) as a paler ring symbol, despite 
uncertainties about is existence. The data points are, with in boldface those 
investigated more closely below: (1) GG Tau Bb, (2) TWA 5 B, (3) GJ 
417 BC, (4) GSC 8047 B/b, (5) DH Tau B/b, (6) GQ Lup b, (7) 2M1207 b, 
(8) AB Pic B/b, (9) LP 261-75 B/b, (10) HD 203030 B/b, (1 1) HN Peg B/b, 
(12) CT Cha b, (13, 14, IS, 16) HR 8799 bcde, (17) Wolf 940 B/b, (18) G 
196-3 B/b, (19) /} Pic b, (20) RXJ1609 B/b, (21) PZ Tel B/b, (22) Ross 
458 C, (23) GSC 6214 B/b, (24) CD-35 2722 B/b, (25) HIP 78530 B/b, 
(26) SR 12 C, (27) HR 7329 B/b, (28) Fomalhaut b. 



and leads to a greater distance between the hot-start cooling curves 
for 10 and 15 Mj than for 5 an d 10 Mj. (This is clearly visible 
in fig. 1 of iBurrows et al.l ( I2001 *). which also shows that there is a 
similar gap for low-mass stars at 10"* Lq and 1-10 Gyr, due to the 
hydrogen main sequence.) In particular, cooling tracks for objects 
of 15, 20, and 25 Mj nearly overlap at ~ 100 Myr around 10"* Lq, 
where data points and their errorbars collect too. This tentative in- 
dication of an agreement between the detections and the cooling 
tracks suggests that the latter might be consistent with the datfl It 
will therefore be interesting to see how significant the gap is and 
how it evolves as data points are added to this diagram. 

We finally note that, in Fig. [8] the best-fitting age is cal- 
culated as the geometric mean of the reported upper and lower 



Cold-start curves too show this gap, to the extent that the luminosity lise 
due to deuteiium burning is very sensitive to the initial entropy (see Fig.|8], 
which would need to be set accordingly finely to have the lightcurves pass 
through the data gap. The implicit assumptions here are that the observed 
distribution of masses is uniform in the approximate range 5-25 Mj (as 
are the mass values chosen for Fig. [8), and that the same applies to the 
initial entropy. The former cannot currently be validated but the latter seems 
reasonable, as the entropy interval over which cooling curves change from 
going above to below the gap is very naiTow. 



bounds since, in most cases, no best-fitting value is provided, and 
the bounds are typically estimates from different methods, which 
cannot be easily combined. In fact, the ages of young (< 5() ( ) Myr) 
stars are in general a challenge to determine, as ISoderblomI (|20I(]|) 
reviews, and represent the main u ncertainty in direct observations. 
Moreover, as lFortnev et al ]^2005) point out, assuming co-evality of 
the companion and its primary may be problematic for the youngest 
objects. Indeed, a formation timescale of ^ I-IO Myr in the core- 
accretion scenario would mean that some data points of Fig.[8]be- 
low ~ 10 Myr may need to be shifted to lower ages, by an un- 
known amount. This consideration is thus particularly relevant for 
GG Tau Bb, DH Tau B/b, GQ Lup b, and CT Cha b (data points 1, 
5, 6, and 12), which are all possibly younger than 1 Myr, and would 
require a closer investigation. 

We now provide detailed constraints for three planetary 
systems, chosen for the low hot-start mass of the companion 
(2MI207 b) or because additional mass information is available 
(HR 8799 and j8 Pic). 



5.2 2M1207 

The companion to the brown dwarf 2MASSWJJ207334-393254 
(2MI207 A, also known as TWA 27 A; Icizisl I2OO2I) is the 
first directly-imaged objec t with a hot-start planetary mass 
jChauvin et ai] |2004 I2OO5I) . Since the age and luminosity of 
2M107 b are the inputs for our analysis, they are discussed in 
some detail in Section lAT] along with tentative information on the 
mass . We adopt an age of 8^3 Myr ( Chauvin et al. 2004; Song et alJ 
2006h and a luminosity of log L/Lq = -4.68 ± 0.05 ( Barman et alj 
201 l b|), and assume that deuterium- burning masses above ^ 13 Mj 
( Spie gel. Burrows. & Milsomll20T ih are excluded. 

Fig. |9] shows the joint constraints on the mass and initial en- 
tropy of 2M1207 b based on its luminosity and ag e. We recover 
the hot-start mass of 3-5 Mj dBatman et alJl20l"Tbh but also find 
solutions at higher masses. If deuterium-burning masses can be 
excluded, the formation of 2M1207 b must have led to an initial 
entropy of 5; ^ 9.2, with an approximate formal uncertainty on 
this lower bound of 0.04 A^s/baryon (see Section [3Tt due solely 
to the luminosity's statistical er ror, in dependent of the age's. This 
initial entropy implies that the cold starts are too cold by 
0.7 /cB/baryon, roughly independently of the mass, to explain the 
formation of t his planet. This is consistent with the time sc ale-based 
conclusion of iLodato, Delgado-Donate. & Clark3 ( l2005i) that core 
accretion catinot be responsible for the formation of this system 
if one also accepts the received wisdom that core accretion neces- 
sarily leads to the coldest starts. However, our robust quantitative 
finding is more general, in that it provide constraints on the result 
of the formation process which are model-independent. 

To show how these M{S \) constraints can easily be made even 
more quantitative and thus suitable for statistical analyses, we ran 
a Metropolis-H astings Markov-Chain Monte Carlo (MCMC; e.g. 
lGregorvll2005l) in mass and entropy with constant priors on these 
quantities. Uncertainties in the luminosity L and age t were in- 
cluded in the calculation of the^^ by randomly choosing an 'ob- 
served' Lobs and a stopping time for the cooling curve fstop at every 
step in the chain. The quantity log Lobs was drawn from a GauBian 
defined by the reported best value and its errorbars, and ?stop from 
a distribution which is either constant in t between the adopted up- 
per and lower limits t\ = 5 Myr and /2 = 12 Myr and zero oth- 
erwise, or lognormal in /, centered at Iq = y/1/2 = 7.7 Myr and 
with (T|og, = Xogtolti = 0.19 dex. The results are shown along the 
vertical axis of Fig.|9]for four different assumptions. To obtain the 
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Figure 9. Allowed values of planet mass M and initial entropy Si 
for 2M1207 b for as sumed ages of 81*^3 Myr The luminosity is that of 



iBarmanetailtollbl) , logLboi/LQ = —4.68 ± 0.05, with the hatched re- 
gions corresponding to curves within ± 1 cr of the luminosity values and the 
thick lines to the central value. The curves along the vertical axis show 
posterior distributions from MCMCs with different mass priors and as- 
sumed luminosity distributions. The solid curves only have an upper mass 
limit of 13 M], while the dashed curves include solutions in the deuterium- 
burning region, up to 20 Mj. For the blue curve (further from the left axis) 
and the corresponding dashed grey posterior (practically identical), a log- 
normal error distribution in time centered at V5 X 12 ^ 7.7 Myr and of 
width 0.19 dex was used, while the red curve and corresponding dashed 
line are for a flat distribution in time betwe en 5 and 1 2 Myr (and zero oth- 
erwise). The circles show the results from Marlev e t al.l ll2007ll . increased 
by 0.38 <:B/baryon (cf. Fig.|4), for cold starts (filled circles) and hot starts 
(open circles), where the entropy is 1 Myr after the onset of cooling. 



two solid lines, we applied an upper mass cut at 13 Mj and took a 
lognormal (less peaked curve) or a top-hat (more peaked) distribu- 
tion in time. The dashed curves come from the same MCMC chains 
but with a mass cut-off of 20 Mj, i.e. including deuterium-burning 
objects. Because of the delayed cooUng due to deuterium burning, 
the required initial entropy drops down faster with mass than in the 
cold-start branch, such that lower 5i values are possible. At 20 Mj, 
the required 5i is 8.2, but it is still only 9 k^lhm-yon at 15 Mj. How- 
ever, the phase space for high masses is very small since the initial 
entropy needs to be extremely finely tuned; hence the smallness of 
the effect on the posteriors. As Fig.|9]shows, different assumptions 
on the luminosity, age, and mass priors all lead to similar results 
for the initial entropy, namely that 5i ^ 9.2 and that there are more 
solutions near this lower limit. 



5.3 HR 8799 

We now turn to the only directly-imaged system with multiple ob- 
jects for which planetary masses are possible, HR 8799. The age of 
the system and the luminosities of the companions are discussed in 
Section lA2l along with information o n the mass. We cons ider ages 
of 20 to 160 Myr, close to the values of lMarois et alj ( l2008l) . and use 
the standard luminosities of logLboi/L p = -5.1 ±0.1 (HR 8799 b ), 
-4.7 ± 0.1 (cd) and -4.7 ± 0.2 (e) from lMarois erai] ( l2008ll2oTol) . 
It also seems resonable to assume that deuterium-burning masses 
can be excluded for all objects, thanks to the (preliminary) results 
from simulations of the system's dynamical stability. 

Fig. [To] shows the joint constraints on the masses and initial 
entropies of HR 8799 b, d, and e. Uncertainties in the age are 
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Figure 10. Allowed values of planet mass M and initial entropy Si for 
HR 8799 b (logLboi/L© = -5.1 ±0.1; blue bands), HR 8799 cd (-4.7 ±0.1; 
purple), and HR 8799e (-4.7 ± 0.2; orange) for assumed ages of 20 (lower- 
mass group) and 160 Myr (higher-mass band; planets cde are outside 
the plotting range). Luminosities are from Maro is et al.l ( I200 8. 2010), and 
nearly in agreement with Marley et al. (2012). The lines correspond to the 
best-fitting luminosity values, while the hatched regions use the l-o" error- 
bars. Along the left axis are shown the posterior distributions in Si (colours 
as for the M-S\ hashed bands), with no particular relative normalisation. 
These were obtained from an MCMC, taking luminosity and age uncertain- 
ties into account based on lognormal distributions (see text), and using flat 
priors in S i and M. The mass upper cut-off is Mn,ax = 1 3 Mj (full lines) or 
'Wmax = (7, 10, 10, 10) Mj (for planets bcde, respectively; dashed lines), 
since analyses of the system's dynamical stability seem to indicate that 
higher masses are unstable. This is indicated by the vertical line segments 
with arrows in the direction of the excluded masses. The circles show the 
results from M07, increased by 0.38 Are/baryon (cf. Fig. |4), for cold starts 
(filled circles) and hot starts (open circles), where the entropy is 1 Myr after 
the onset of cooling. 



taken into account by considering the two extremes of 20 Myr and 
160 Myr separately, while the 1-cr errors in the luminosities are re- 
flected by the width of the hatched regions. We find hot-start masses 
for 20 Myr of 4.4^°;^ (b), 6.3 ± 0.6 (cd), and 6 ± 1 Mj (e), where 
the errorbars here come only from t hose on the luminosity . These 
values are in good agreement with iNeuhiiuser & Schmidtl Jioilh 
and are sim ilar to the usually-cited 30-Myr values of (5,7, 7, 7) Mj 
jMarois et al. 2010) . The hot-start masses for 160 Myr are above 12 
(HR 8799 b) and 13 Mj (cde), and the M{S{) constraints for the lat- 
ter three are not shown since they are within the deuterium-burning 
regime. 

Excluding deuterium-burning masses for all objects and using 
only the luminosity measurements, Fig. [10] shows that all planets 
of the HR 8799 system must have formed with an initial entropy 
greater than 9 A;B/baryon, with 5i > 8.9 for b, 5i > 9.2 for c 
and d, and 5i > 9.1 for e. (When comparing models of different 
groups, one should remember that there is typically a systematic 
uncertainty in entropy values of A5syst - 0.4, as mentioned in Sec- 



tion |2.4| ) Using tentative upper mass limits of 7 and 10 Mj, respec- 
tively, the lower bounds on the initial entropies can be raised to 9.2 
(b) and 9.3 A-B/baryon (cde) if one takes the conservative scenario 
of the l-cr lower luminosities value at 20 Myr These lower bounds 
on the entropy are however mostly independent of the age because 
they are set by cold-start solutions, where the age is much smaller 
than the cooling time at that entropy. Formal uncertainties on the 
lower bounds due to those in the luminosities are (see Section lSTt 
approximately 0.07 (bed) or 0. 14 ^b /baryon (e) and thus negligible. 
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Here too we ran an MCMC to derive quantitative constraints 
on the initial entropy of each planet. The quantity log Lobs was 
drawn from a GauBian defined by the reported best value and its 
errorbars, and fstop from a distribution which is either constant in 
t between the adopted upper and lower limits ti = 20 Myr and 
t2 = 160 Myr and zero otherwise, or lognormal in t, centered at 
'o = '^hh = 57 Myr and with crioe, = logf2/fi = 0.9 dex. Pos- 
teriors on the initial entropy for each of the HR 8799 planets are 
shown along the vertical axis in Fig. (TO] using a flat prior in 5i 
and a mass prior constant up to an Mn,ax and zero afterward. The 
cases 'without mass constraints' (Mmux = 13 Mj) are nearly con- 
stant in Si, especially for planet b, but show a peak near cold- 
branch values of 9 and 9.5 ^re/baryon for HR 8799 b and cde. 
Adding mass information from dynamical-stability simulations by 
taking M„,„ = (7, 10, 10, 10) Mj flattens the S i posterior and shifts 
the minimum bounds at half-maximum from (8.9, 9.2, 9.2, 9.2) to 
(9.3,9.5,9.5,9.4) A;B/baryon, respectively. We note that these re- 
sults are insensitive to both the form of the uncertainty in t and 
the use of a non-flat prior in mass (as shown below for /3 Pic b in 
SectionlSH. 

Comparing to the 'tuning fork' entropy value s repr oduced 
in Fig. [Tol we find that the coldest-start models of IM07| cannot 
explain the lumino s ity m easur ements for the HR 8799 planets. 
ISpiegel & Burrowsl ( 1201 2|) and iMarlev et aO ( 1201 2h also came to 
this conclusion, with the latter noting that 'warm starts' match 
the lumi nosity constr a ints. I t is now possible to say specifically 
that the iMarlev et al.\ i200'^} cold starts would need to be made 
AS - 0.5 kglbaryon hotter to explain the formation of the HR 8799 
planets. Given that the precise outcomes of the core accretion and 
gravitational instability scenarios are unc ertain and that this sys- 
tem represents a challenge for both (as iMarois et al.l |2010| and 
ICurrie et al.ll201 ll review), quantitative comparisons such as our 
procedure allows should be welcome to help evaluate the plausi- 
bility of the one or the other 

We note in passing that one needs to tak e care also 
when interpreti n g the measurements of iHinklev et al.l (1201 lah and 
IClose & Malel j2010 l). These authors measured upper limits on 
the brightness of companions within 10 AU and between 200 and 
600 AU from the star, resp ectively. How e ver, b oth groups then 
used the hot-start models of iBaraffe et al.l j2003h to translate the 
brightness limits into masses (11 Mj at 3-10 AU and 3 Mj within 
600 AU). Therefore, since colder-start companions would need to 
be more massive to have the same luminosity, what they provide are 
really lower upper limits on the mass of possible companions. How 
much higher the masses could realistically be in this case is hard to 
estimate without a bolometric luminosity, but there is an important 
general point: without the restriction of considering only hot-start 
evolutionary tracks, luminosity upper limits do not provide unam- 
biguous mass constraints. Incidenta lly, this more general view of 
the results of iHinklev et alj J2OI lal) means that the unseen com- 
panion evoked bv lsii et alj j2009| ) as the possible cause of the inner 
hole (at < 6 AU) does not have to be of small mass. This inner 
object would nevertheless have to be consistent with the results of 
dynamical stability simulations. 
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Figure 11. Allowed values of planet mass M and initial entropy Si for 



/? Pic b using our estimated logLboi/L© 



-3.90;""5 and an age of 



12+= IVIyr. The hatched regions correspond to luminosity values within the 
(asymmetrical) 1-cr interval, with the central values marked by thick lines. 
The curves along the vertical show the marginalised posterior distribution 
on Si from an MCMC simulation using the luminosity and age values and 
their uncertainties. These were taken to be lognormal but asymmetric for the 
luminosity, and lognormal in time, centered at ?o = V8 X 20 = 12.7 Myr 
with a standai'd deviation of 0.35 dex. The full blue line (closer to the 
vertical axis) results from the (M, Si) distribution being multiplied by a 
dNJdM oc M" prior, with o: = - 1 .3 iCumming et al.l2008| - |Nielsen & Clos3 
I2OIOI) . and zero for M > 12 Mj, while the full grey posterior has only a mass 
cut at 12 Mj, from the iLagrange et all l |201 2b) radial- velocity constraint. 
The dashed and dotted curves both have an upper mass cut of 20 Mj , with 
and without the d N/dM prior applied, respectively. The circles show the 
results from Marl ev et al.l l l2007 1. increased by 0.38 /ia/baryon (cf. Fig.|4j, 
for cold starts (filled circles) and hot starts (open circles), where the entropy 
is 1 Myr after the onset of cooling. 



detail in Section |A3l how we derive a bolometric luminositj^ of 
logLboi/Lo = -3.90:«;»^. 

Fig.ll llshows the M(Si) constraints available for /? Pic b from 
our luminosity estimate and the radial-velocity (RV) constraint. 
We recover a hot- s tart mass ^ 9.5 ± 2.5 Mj, i n agr eement with 
lOuanz et all ^2010) and Neuhauser & Schmidt! (^2), but addi- 
tionally find that higher masses are consistent with the luminosity 
measurement. Excluding solutions where deuterium burning plays 
an important role in the evolution of the object (recognisable by the 
extreme thinness of the constant-luminosity M(Si) curve) implies 
that Si ^ 9.8. Using the RV mass upper limit, these constraints 
on the initial entropy can be made tighter: with an age of 12 Myr, 
it must be that Si J? 10.5. Since this corresponds to a warm start, 
both uncertainties on the age and on the luminosity contribute to 
that on the minimum Si, of order 0.5 /^B/baryon. 

This lower limit on Si implies that coldest-start objects of any 
mass are too cold by a significa nt 1.5-2.0 ks/hary on. Various au- 
thors (e.g. lOuanzet al . 2010; B onnefov et alJ201l|k recognised that 
the classical cold starts ( M07I : IFortnev et al.ll2008l) cannot explain 
the observations, but it now possible to quantify this. These results 



5.4 J3 Pic 

A co mpanion to the well-studied star P Pic was fi rst observed in 
2009 dLagrange et"ZI 120091; Ifionnefov et aH I2OI ih and, very re- 
cently, it became the first directly-detected object with a planetary 
mass for which radial- velocity data is also av ailable. The age of 
the system is 12;^* Myr dZuckerman et aklboOlh . and we discuss in 



We note that as this manuscript was being prepared, w e became aware 
of the first robust estimate of the bolometric luminosity, by Bonne fov et al] 
(2013). They find logLboi/LQ = -3.87 ± 0.08, which excellently agrees 
with our value and thus does not change our conclusions. In particular, they 
find similai' constraints on the initial entropy of /? Pic b, although this de- 
pends on which band they use (cf. their fig. 1 1 with our Fig. ll 11 . 
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are mostly insensitive to the uncertainty on the age range; using in- 
stead 12-22 Myr as summarised by Femandez^t al. (2008) would 
not change the conclusions. 

As for the objects in the 2M1207 and HR 8799 systems, we 
ran an MCMC to obtain a posterior distribution on the initial en- 
tropy. This is shown along the vertical axis of Fig. [TT] for four 
different assumptions. In all cases, we assumed lognormal uncer- 
tainties on the age and the luminosity, with asymmetric upper and 
lower errorbars for the latter For the full curve closer to the verti- 
cal axis (in blue), we applied an upper mass cut at 12 Mj and took 
into account that the underlying (real) mass distribution is possibly 
bias ed towards lower masses, as radial velocity me asurements indi- 
cate jCumming et alj2008l : lNielsen & Closjioid) . Out of simplic- 
ity, this was done by taking the (M, Si) distribution obtained with a 
flat prior in mass and multiplying it by a AN I AM oc M" prior, with 
a = -1.3, thus weighing lower masses more. Of course, the value 
of a might depend on the formation mechanism relevant to the ob- 
ject but this serves to illustrate the effects of a non-constant prior 
on mass. The other three posterior distributions on 5i of Fig.fTTIdn 
grey) correspond to the remaining combinations of 'with mass cut 
or not' and 'with power-law mass prior or not' . These curves are all 
similar, with the radial-velocity measurement increasing the mini- 
mum bound at half maximum from 10.2 to 10.5 ^a/baryon, quite 
insensitively to the use of the AN I AM prior 

At the distance from its primary where p Pic b is cur- 
rently located (- 9 AU), core accretion is expec ted to be effi- 
cient and thus a likely mecha nism for its formation ( Lagrange et alj 
1 201 ll : iBonnefov et alj[2013h . Thus, the question posed to forma- 
tion models is whether core accretion can be made hotter (by 1 .5- 
2 /^e/baryon) than what traditional cold starts predict. Since these 
coldest starts assume that all the accretion energy is radiated away 
at the shock, the constraints on the initial entropy stress the need to 
investigate the physics of the shock, when the initial energy content 
of the planet is claimed to be set. 



6 SUMMARY 

The entropy of a gas giant planet immediately following its for- 
mation is a key p arameter that can b e used to distinguish planet 
formation models jMarlev et al .1120071) . In this paper, we have ex- 
plored the constraints on the initial entropy that can be obtained for 
directly-detected exoplanets with a measured bolometric luminos- 
ity and age. When the initial entropy is assumed to be very large, a 
'hot start' evolution, the measured luminosity and age translate into 
a constraint on the planet mass. In contrast, when a range of initial 
entropies are considered ('cold starts' or 'warm starts'), the hot- 
start mass is in fact only a lower limit on the planet mass: larger- 
mass planets with lower initial entropies can also reproduce a given 
observed luminosity and age. Fig. |5] shows the allowed values of 
mass and entropy for different ages and luminosities, and can be 
used to quickly obtain estimates of mass and initial entropy for any 
given system. 

To derive these constraints, we constructed a grid of gas gi- 
ant models as a function of mass and central entropy which can 
then be stepped through to calculate the time evolution of a given 
planet. In a hot-start evolution, the structure and luminosity of cool- 
ing gas giant planets is usually thought of as being a function of 
mass and time only. Once the assumption of hot initial conditions 
is removed, however, a more convenient variable is the entropy of 
the planet. One way to think of this is that gas giants obey a Vogt- 
Russell theorem in which the internal structure, luminosity, and ra- 



dius of a planet depend only on its mass and entropy. Fig.|2]shows 
the luminosity as a function of entropy for different masses, and a 
general fitting formula for L(M,S) is given in Section |T2l A note- 
worthy result is that in the intermediate-entropy regime (5 - 8.5- 
10 A:b /baryon), where the outer radiative zone is thick and follows a 
radiative-zero solution, the luminosity obeys L = M f{S ) as found 
for irradiated gas giants by lArras & Bildstenl (|2006|) , with f(S ) a 
steeply increasing function of the entropy. 

We find that our models are in good agreement (with i n tens 
of p ercent) with t he ho t-start models of iBurrows et al.l ( Il997h 
and iBaraffe et aL N2003h . as well as the cold-start models of 
iMarlev et alj ( l2007h . a nd cooling models c alculated with the mesa 
stellar evolution code dPaxton et al.ll20lil 12013). However, there 
is a small offset in the entropy of =^ 0.4 /baryon between the 
m odels calculated here a n d using mesa o n the one h and and those 
of iBurrows elal] d 1 9971) . iBarafi'e et al.1 j2003l) , and iMarlev et alj 
( l2007l) on the ot her. Thi s may be due to a difference in the imple- 
mentation of the lSCvril equation of state. Uncertainties due to dif- 
ferences in opacities or composition do not appear to be en ough to 
account for the offset; for example, we estimate from .Saumon et al.l 
( Il996l) that the uncertainty in the helium (Y) and metal (Z) mass 
fractions introduces variations of at most 10 per cent in the lumi- 
nosity. 

We stress again that when the initial entropy is allowed to take 
a range of values, the hot-start mass is only a lower limit on the 
companion mass. The larger range of allowed masses means that 
the hot-start mass could actually lead to the mischaracterization of 
an object, with a hot-start mass in the planetary regime actually cor- 
responding to an object with a mass above the deuterium-burning 
limit for low enough entropies. One way to break the degeneracy 
between mass and entropy is an accurate determination of the ra- 
dius from spectral fitting (or actually a determination of log^ and 
T^ff from the spectrum), which would yield the mass and (current) 
entropy of the object without any degeneracy (see e.g. Fig. |7](. As 
discussed in Section O however, current atmosphere models have 
significant uncertainties that make this approach difficult. Another 
possibility is to obtain independent constraints on the mass of a 
companion, for example from dynamical considerations. 

In Section |5] we applied our models to three directly-imaged 
objects which have hot-start masses in the planetary-mass regime 
and for some of which there are additional constraints on the planet 
mass. We find that the initial entropy of 2M1207 b is at least 
9.2 /baryon, assuming that it does not bum deuterium. For the 
planets of the HR 8799 system, we infer that they must have formed 
with 5i > 9.2 /tb /baryon, independent of the age uncertainties for 
the star Finally, a similar analysis for /? Pic b reveals that it must 
have formed with 5i > 10.5 /baryon, using the radial-velocity 
mass upper limit of 12 Mj. These initial entropy values are respec- 
tively ca. 0.7, 0.5, and 1.5 fcg /b aryon higher t han th e ones obtained 
from core accretion models bv IMarlev et al.l (l2007b . This quantita- 
tively rules out the coldest starts for these objects and constrains 
warm starts, especially for /3 Pic b. 

An important point is that the uncertainties in age and lumi- 
nosity impact the derived hot-start mass and the lower bound on 
initial entropy in different ways. The major uncertainty in direct 
detections is the age of the star and, relevant for very young sys- 
tems, all the more that of the planet. This age uncertainty translates 
into errorbars for the hot-start mass which are AM/M ^ \At/t. 
However, the initial entropy on the mass-independent branch of the 
M(Si) curve is due only to that in the bolometric luminosity, with 
A5 i - 1 /AA log Lboi where A =^ 0.7 when 5 i < 9.5. This uncertainty 
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AS j is typically very small, which means that, up to systematic er- 
rors, the initial entropy can be determined quite accurately. 

It has been pointed out before that current direct-imaging de- 
tections are all inconsistent with the cold-start predictions from 
core-accretion models. It is important to stress however that the 
cold starts are in some sense an extreme case, as they assume 
complete radiation efficiency at the shock during runaway accre- 
tion, which is argued to set the low initial entropy of planet. By 
varyin g the nebula temperature or the accretion rate, iMarlev et alj 
( I2OO7I) were able to chang e the entropy by barely AS - 0.1, 
but Mordasini et alj ( 1201 2bl) report that it is possible within the 
core-accretion scenario to obtain considerably higher entropies, in- 
creased by as much as 1-2 fcB/baryon (Mordasini et al., in prep.). 
With y6 Pic b a likely candidate for formation by core accretion, this 
indicates that it is essential to gain a deeper understanding about 
what sets the initial entropy, for instance by looking in more detail 
at the properties of the shock during runaway accretion. 

The derived bounds on S \ for the HR 8799 objects and p Pic b 
made use of information on the mass, which comes from dynamical 
stability analyses and radial velocity, respectively. Radial-velocity 
data of directly-detected planets is currently available only for 
jS Pic b, but this should change in a near future as close-in planets 
start being detected directly. In the absence of dynamical informa- 
tion, an upper limit to the mass (and thus a lower limit on the initial 
entropy) should be obtainable from log g, even if its errorbars are 
large; thus, in practice, the M{S\) constraint curve does not extend 
to arbitrarily high masses as a pure luminosity measurement would 
imply. 

Finally, we ran Markov Chain Monte Carlo simulations to de- 
rive more detailed quantitative constraints on the mass and entropy 
of directly-detected objects. When taking the uncertainties in the 
age and luminosity into account, we considered normal, lognor- 
mal, and flat distributions and found the chosen form to make little 
difference. The advantage of this approach is that it allows one to 
derive posterior distributions on the initial entropy, which are suit- 
able for statistical comparisons. Given the small semi-major axis 
of yS Pic b (9 AU), we also tried a prior dN/dM oc M^'-^ (which 
describes the population of radial-velocity planets) in addition to a 
flat prior on M. The latter case yielded posterior distributions on Si 
with a more pronounced peak. 

The benefits that the expected large increase in the number 
of directly-detected exoplanets in the near future should bring are 
at least twofold. Firstly, each new detection will yield a new con- 
straint on the initial entropy and therefore formation mechanism. 
Particularly with the ability to detect lower mass gas giants at small 
semi-major axes (< 20 AU), which should be possible with in- 
struments such as GPI or SPHERE, there will be an opportunity to 
constrain the state of the gas giant immediately following core ac- 
cretion. Secondly, with a larger sample of objects comes a chance 
for a statistical comparison with formation models. For example, 
as noted in Section [STI the onset of deuterium burning leads to a 
relative underdensity of planets with luminosities ~ 10"* Lq and 
ages of tens of Myrs (assuming that the masses of substellar objects 
are smoothly distributed near the deuterium-burning limit). Indeed, 
it is interesting that there appears to be such an underdensity in 
the current data sample (see Fig. 8), although the small number of 
detections so far means that this could be due to a statistical fluctu- 
ation. The best constraints on initial conditions for planet cooling 
will come from improved spectral models that can give reliable de- 
terminations of log g and T^g. 
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APPENDIX A: AGE, LUMINOSITY, AND MASS 
CONSTRAINTS 

Al 2M1207 

Al.l Age and luminosity 

The first directly-imaged object with a hot-start planetary mass 
JChauvin et aP |2004 l2005h i s located 8" from its M8 brown-dwaif 
primary, a well-stu died (see iGizisI |2002| for the discovery paper and 
ISkemer et al.l l201ll for a summary) membe r of the young ( 8^j Myr) 
TW Hydra; association iChauvin et al.ll2004l: ISong et al.ll200db at a dis- 
tance of 53 ± 1 pc (as avera ged by Skemer et al. 2011 from Gizis 2002; 
Gizis et al.l2007tlBiller & Close.2007.; Mamaiek.2005. : Mamaiek & Meyer. 



2007t iDucourant et al.l |2008|) . implying a projected orbital separation of 
41 ± 1 AU. 

Determining the lumino sity of 2M1207 b is no t straightforward. Pho- 
tometry is available in the 7 iMohantv et al.l \2 00i^. H, K, , and L' bands 
from NaCo observations at the VLT l lChauvin et al. 2004; iMohantv et alj 
l2007h . at 0.8-1.0, 1.0-1.2, 1.35- 1.55, and 1.4-1. 8 |xm (similar to H) 
from the Hubble Space Telesco pe jSong et alj|2 00dl). an d in the Hersche l 
SPIRE bands of 250 and 350 um lRiaz et alj J2012b a). Ske mer et aklfeOllh 
also obtained an upper limit at 8.3-9.1 yim. Spe c troscopy is available 
at 1.1-1.35 Hm, 1.4 -1.8 \xm jChauvin et ai]l2004 IMohantv et alj|2007l ; 
iPatience etal]|2010l). and 1 .95-2.5 |xm JPatience et alj |2010|) and 2.0- 
2.4 um jMohantv et al.l2007l) . i.e. in JHK g and slightly redward . Summing 
up the fluxes in /, H, K,, and L' listed in lMohantv et alj J2007h and using 



the known distance gives a luminosity of 2.1 X 10^^ erg s"', or 27-30 per 
cent of a bolometric luminosity log Ltoi/LQ = -4.68 to -4.74 (see below). 
A 'direct luminosity dete rmination' is thus not possible. 

As summaiised by iBarman et al Itj) , there are inconsistencies 
between the luminos ity of 2M1207 b and its spectral type, determined by 
IChauvin et all J2004 

to be L5-L9.5. With the corresponding bolometric 
correction BCy of loolimowski et al ] J2004 . t he Kg-banA magnitud e im- 
plies a luminosity of log L/L© = -4.7 ± 0.1 ' Barman et ^1201 Ihl) . The 
hot-start, equilibrium models of iBarafFe et alj ■ 20031) then yield from the 
age a nd luminosity an effective temperature of 1010 ± 80 K IBarman et alj 
[201 lb) in disagreement with Tefi- 1600 K implied by the spectral type 
(Mohanty et al. 2007; Patience et al. 2010, 2012). However, this second 
value is questionable. Indeed, the high effective temperature and low lu- 
minosity would require an improbably small radius of =: 0.6 Rj. 

There are two distinct approaches to the solution. IMohantv et alj 

i2007i) argue that the actual luminosity has been underestimated due to 
grey extinction by an al most edge-on disc (the ^ 25-Mj brown -dwarf pri- 
mary is accreting; e.g. IStelzer, Scholz, & Java wardhana '2007'). The cor- 
respondingly higher luminosity is reported by Ducourant et al. (2008) as 
logLboi/L© = -3.8 ± 0.1. On the contrary, Ske mer et alj ( 1201 1) argue 
against the disc explanation and sug gest th at thicker clouds than what might 
be naively expectej*^ are needed. IBarman et alj l j201 lalbl) go further and 
explicitly claim that the problem is with the derived T^g. They show that 
an atmosphere model with T^g = 1000 K can fit very well the photome- 
try and spectroscopy if clouds of typical thickness and also, crucially, non- 
equilibrium chemisty are included. The latter leads to a heavily reduced 
methane abundance (by ca. two orders of magnitude at photospheric depths) 
compared to the chemical-equilibrium clouds and thus to redder co lours 
than expected. From their best-fitting models, IBarman et alj 1201 Ihl) esti- 
mate a luminosity of log L/Lq = -4.68 ± 0.05 , in agreement with the lu- 
minosity derived from BC^: . As lLuhmMl j2012h notes, differing bolometric 
corrections for old (field) and young brown dwarfs are thus not the sole 

explanation. 

We adopt the luminosity of lBarman et alj J201 Ibl) , in agreement with 
INeuhauser & Schmidll J2OI2I) who report log L/L© = -4.74 ± 0.06 (with 
however Teff = 1 590 ± 280 K). Using the age of 5-12 Myr and the 
iBaraffe et alj J2003l) cooling t r acks, th is lumi nosity yields a hot-st art mass 
of 2-5 M] as ' Barman et alj j2011hl) state. I Skemer et all fcOllI) give a 
slightly higher hot-start mass between 5 and 7 Mj based on T^g = 1000 K 



See comment in section 1.2.2 of lMarlev et alj fe012l) . 



and the iBurrows et al ] 1 I1997I) models. We note that these errorbars match 
the estimate AM/M ^ ^^t|t ^ 0.5 from SectionlJT] 



Al.l Mass information 

Mass information of dynamical origin for this two-body system is not avail- 
able. Indeed, since the separation implies a period of at least (depending on 
the eccentricity) 1700 yr, detectable orbital evolutio n or change of th e ve- 
locity amplitude are not expected in the near future j Mamaie3l2005h . pre- 
cluding both astrometry and radial-velocity measurements. 

However, the surface gravity is somewhat cons trained, which c an be 
used to set an approximate mass upper limit. Barma net al.l ( l2011bh state 
that their best-fitting model has g = Iff* cms"^ but do not provide any 
sense of how large the uncertainty on this value migh t be. However, typica l 
errorbars (as in then' si milar analysis for HR 8799 b; IBarman et al.ll201 lal) 
are at least of 0.5 dex. IMohantv et al ] J2007l) found that the fit to both the 
photometric and spectroscopic data is rather insensitive to log g within 3.5- 
4.5 (cms^^), in agreement with the indications of low gravity from the tri- 
angular //-band spectral shape and rel atively weak Na 1 absorption (e.g. 
lAllers et alj|2007l; IMohantv et alj |200"7l) . Also. 'Patien ce et alj iioH) fit J, 
H, and K spectra w ith five grids of atmosphere mode ls, including BT -Settl 
jAllard, Homeier, & Freytag 2011), Drift-PHOEmx <Helhng et akbOO S*). and 
those of Marley et al. ( Ackennan & Marley 2001). The Marley et al. models 
yielded a gravity on the edge of their grid {g = 10^" cms"^), but the others 
gave logg ^ 3.5, 3.5, 4.3, and 5.0 (cms^^), respectively. (With the best- 
fitting Teff = 1500-1650 K, the implied radii are of 0.4-0.7 Rj, well below 
any theoretical cooling track.) Even though there are systematic issues with 
atmosphere models of young, low-mass brown dwarf's are expected, we will 
take these results to suggest tentatively that 2M 1207 b has a low gravity. 
As discussed below, the initial entropy on the cold-start branch is Si = 9.2, 
which is thus an upper limit to the current entropy. With S < 9.2 and a 
reasonable upper limit logg < 4.0 (cms^^), the upper bound on the mass 
is ^ 7 Mj, and for logg = 4.35 (cms^^) it is 12.7 Mj. Therefore, we shall 
assume that the estimates of the surface grav ity imply a mass b elow the 
deuterium-burning limit, which is near 13 Mj jSpiegel etal.l201ll) . 



A2 HR 8799 

A2. 1 Age and luminosities 

Several properties of HR 8799 let its age be estimated: variability from 
non-radial oscillations, low abundan ce of iron-peak el ements, and far- 
IR excess due to crrcumstellar dust jMarois et alj|2008l) . Along with its 
Galactic space motion and p osition in a Hertzsprung-Russell diagram, 
these let Marois Tt alj J2008l) estimate an age range of 30-160 Myr 
with a prefeiTed value of 60 Myr, consistent with the 20-150 Myr 
range of Mooret al. (2006) based on membership in the Local Associa- 
tion. Recently, .Baines et aL 120120 used interferometric measurements of 
HR 8799's radius to derive a stellar mass and age. They found best-fitting 
ages of 33 or 90 Myr, depending on whether the star is approaching or mov- 
ing away from the main sequence. However, the statistical errorbars, which 
do not take uncertainties in the stellar models into account, are considerable 
in the second cas e (the l-cr ran ges are 20^0 and 40^71 Myr, respectively). 
Nevertheless, if iBaines et akN measurement of the stellar radius and the 
deduced metallicity are corre ct, the age range of 1.1-1.6 Gyr from the as- 
te roseismologic al analysis of Mova et al] j2010al) would be compromised, 
as IBaines et al] 12012) point out. Indeed, they estimate a near- solar metal- 
Ucity, w hich contrasts with the [M/H] =: -0.3 or -0.1 result of lMova et alj 
J2010bh . while metallicity is an important input of asteroseismolog ical anal- 
yses. M oreover, there is the statistical argument put forth by M arois et alj 
Hooi) that massive discs (such as HR 8799's of 0.1 Mm; iSu et aU2009l) are 
unlikely to be found around older stars. In our analysis, we shall therefore 
ignore the 1.1 Gyr result and instead use 20 and 160 Myr as lower and up- 
per lim its, which brackets the ranges reviewed in Moya et al. (2010a) and 

IBaines et al. (2012)^ 

[Marois et ah, J2008l I2OIOI) estimated luminosities of logL(,oi/LQ = 
-5.1 ± 0.1 (HR 8799 b), -4.7 ± 0.1 (cd) and -4.7 ± 0.2 (e) from the known 
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distance of 39.4 ±0.1 pc and six infrared magnitudes, covering ^ 40 per 
cent of the bolometric luminosity, and also^from bolometric corrections for 
brown dwarfs. For their part, iMarlev et al.l (|2012j) recently derived lumi- 
nosities of log Lboi/L© = -4.95 ± 0.06, -4.90 ± 0. 10, and -4.80 ± 0.09 for 
planets b, a*H and d, respectively, by self-consistently obtaining the radius 
from evolutionary models. This contrasts with the usual procedure of opti- 
mising (R/d)^ along with T^g and log g to fit the photometry , which yields 
unphy sically small radii o f ^ 0.8 Rj ^Barman et al.ll201 la ; IMarlev et alj 
I2OI2I) . Since iMaiiev et all considered only hot starts, i.e. fixed 5; = high, 
the cooling tracks are Lba\(M, f) and R(M, t) relations, which let R and t be 
uniquely detemiined from log g and Teff . This gives an age of 360 Myr for 
HR 8799 b and age ranges of 40-100 Myr and 30-100 Myr for c and d, 
consistent with other literature estim ates. Given the difficulties in obtaining 
a reasonable fit. lMarlev et al.l J2OI2I) warn that the first result should not be 
taken seriously, and note that the tension in the age would be reduced by 
c onsidering colder initial co nditions. In our analysis, the luminosity values 
of lMarois et al] )2008ll2oTol) will be u sed since they are sta ndard and almost 

or marginally consistent with those of lMarlev et al.l J2012l). 

Finally, an age of 30 or 60 Myr for the system lead Marois et al. ( 20^ 
to derive from the luminosity and the cooling tracks of Baraft'e et al. ( 20Q. 
masses of (5, 7, 7, 7) or (7, 10, 10, 10) Mj. The uncertainty in the age implies 
(see Section|3) AM ^ 1-1.3 Mj for the hot-start values. 



A2.2 Dynamical stability 

Since it is to date the only directly-imaged multiple-planet system, 
HR 8799 has received a considerable amount of attention with reg ard to 
its dynamical stability (e .g. Marois et al. 2008; Reidemeister e t al]l2009l; 
Fabrvckv & Murrav-Ciavl BOIQ; Moro-Martm et al. 2010,; Marois et alj 



r vckv & Murrav-Ulavl 1 

201oj;ICurrie et al.l201ll;ISudol & H aghighipour 2012 tlEsposito et alj2013F 
Currie et al.l l2012bl; see review m .Sudol & Haghighipouj |2012|) . How- 
ever, only the most recent studies were able to consider all four 
planets. Crucial questions include whether there are two- or three- 
planet mean-motion reson ances (as Gozdziewski & Migaszewski 2009 and 
iFabrvckv & Murrav-Clavll20 10 suggest), what the inclination and eccen- 
tricities of the orbits are (for instance, Lafreniere et al. 2009 estimated 13- 
23° for the inclination of HR 8799 b with res pect to the plane of the sky) 
and whether they are co-planar (against which ICurrie et al 1l2012bl recently 
provided evidence), and, naturally, what the masses (including that of the 
star) are and how long the system should be required to survive. Solutions 
are ver y sensitive to these pa rameters and even to the numerical integrator 
used, as lEsposito et al.l i2013l) note. The parameter space's high dimension- 
ality makes a proper exploration - i.e. without artificially -imposed restric- 
tions as all authors had to assume -, computationally prohibitive, and trying 
to include information about the disc would only make matters worse. 

Stability is estimated by using astrometric constraints and numer- 
ically evolving the system over time, requiring that it be stable (with- 
out collisions nor ejections) for a pe riod equal to its age. However , 
iGozdziewski & Migaszewski ^20091) and lFabrvckv & MuiTav- ci^ feOlOl) 
point out that if it is young with respect to its main-sequence lifetime, 
HR 8799 could indeed be a transient system undergoing dynamical relax- 
ation. Therefore, it may not be possible to draw firm conclusions even from 
the results of a complete analysis. 

Nevertheless, if the direction in which these studies point is coiTect, 
the planets should have as low masses as allowed, with however higher 
masses permitted if some orbits are resonant. For this reason, we shall con- 
sider as approximative upper masses from stability analyses 7, 10, 10, and 
10 Mj for HR 8799 b, c, d, and e, respectively. In particular', while the oppo- 
site cannot be soUdly excluded, it seems plausible that none is a deuterium- 
burning object. 



In fact, the en'orbars on the luminosity of HR 8799 c are clearly non- 
GauBian, but this will not be taken into account out of simplicity. 



A3 ySPic 

A3.1 Age and luminosity 

The namesake A5 dwarf of the nearby (9-73 pc; \Mak) et al 1'2013|) B Pic- 
toris moving group has an age of 12^^ Myr IZuck erman et alj 1200 ll) 
and asymmetric outer and warped inner discs, which have been observed 
for more than tw o decades ( Lagrange et al. 2009, 20 12aL see review in 
Lagrange et a l. 2011 ). A companion was first detected in L' (Lagran ge et alj 



and subsequent ly confirmed a t 4 [xm jOuan z et al. 2010) and in 
I Bonn efov et al.l201ll) . Very recentlv. lBonnefov et aljj2013b added to these 
observations photometry mJ,H, and M' . The distanc e of 19.44 ± 0.05 ptp^ 
to 13 Pic b implies an orbital separation of 8—9 AU iChauvin et al I I2OI2F 
iBonnefov eraUl2013l) . which is the smallest of all low-mass directly- 
detected objects. A n object at this position had been predicted from the 
disc morphology bv lFreistetter. Krivov. & Lohnd j2007l) . 

Until recently, the only bolom ettic luminosity estimate for /3 Pic b is 
due to lNeuhauser & Schmidll J2012t) . who report log Z-boi/Lo = -3.90 



'Y 1 

They fi rstly derived, from the Teff = 1700 ± 300 K value of lBonnefov et alj 
i201ll), a spectral type SpT ^ L2-T4 using the Tj-ff-SpT relation of 
iGolimowski et alj i2004l hereafter IG04|) . They then estimated from their 
SpT-BCs: curve a bolometric correction BCk = 3.3^| qq (T. Schmidt, 
priv. comm.). Thus, the large, asymmetric lower errorbar on the luminosity 
comes from the large, asymmetric lower ertorbar on the bolometric correc- 
ti on, wh ich itself is due to the flat Teff-SpT relation between L7 and T4 
in lG04l A more direct approach to the bolometric luminosity consists of 
converting the colour to a spectral type and obtaining from this a bolomet- 
ric correction. Also using the fits of G0 4, this gives the same luminosity 
as found by Neuhauser & Schmidll j2012h but with a smaller lower errorbar 
of 0.12 dex. One sho uld ho wever note that BCa:(SpT) is not a monotonic 
function (see fig. 6a of lG04h . so that the errorbars are strongly non-Gau6ian. 
With a maximum BCa: of 3.3 mag near L4, log Z^bQl/L n can not formally 
be above -3.9 dex. Howe ver, the BCa:-SpT re lation of IG04| was derived 
for field dwarfs (see also IStephens et al and the spectral classifi- 

cation of you ng objects is not yet well un derstood nor, in fact, well de- 
fined (see e.g. lLiu et al ] |201ll ; lFahert:v et alJi2012 ). Keeping in mind these 
uncertainties in interpreting the photometry, we shall use for the analysis 
log L|,oi/Lq = -3.90^[J"p, where the upper errorbar reflects the residuals of 
the.GQ4fit. 



A3. 2 Mass information 

The object /? Pic b is particularly interesting because it is the first directly- 
imaged companion for which radial-velocity (RV) data is also availablj'**! 
iLagrange et al.ll2"012bl) . Using new and archival data spanni ng eight years 
and t h anks to the high inclin ation of the system (88 ± 2° ; IChauvin et alj 
l2012ll . lLagrange et alji2012bh were able to place tentative lower mass lim- 
its of 1-2 Mj, which is fully consistent with all reasonable age and luminos- 
ity combi nations, even all owing for very large lower errorbars on the latter. 
However. iLagrange et all s u pper limit of 10-25 Mj, with 12 Mj for the 
most probable orbit of 9 AU iLagrange et al. 2009; Ch auvin et al ll2012h IS 
an important result which excludes high-mass solutions and puts the object 
quite likely in the planetary (non-deuterium burning) range. 



*^ This is the value obtained from a re-reduction of Hipparcos data by 
Ivan LeeuwenlJioOTi). However, a n umber of recent studies still use the value 
of 19.3 ± 0.2 PC iCrifo et al.l 19971) . 

To this rarity contribute both intrinsic detection biases - direct imaging 
favours planets further out from their star, resulting in a small RV signal 
- as well as selection biases - target stars are usually chosen based on the 
presence of a disc, which implies that the systems are preferentially seen 
face-on. A possible further hindrance is that young stars - young systems 
being of more interest - are usually active and thus less amenable to radial- 
velocity measurements. 
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